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SUMMARY 


To  establish  basic  definitions  and  to  serve  as  a  reference  for 
future  theoretical  and  applied  research  efforts,  the  fundamental  con¬ 
cepts  of  the  theory  of  continuous  mass  media  are  reviewed,  and  the 
mathematical  foundations  of  constitutive  equations  for  isotropic  and 
homogeneous  materials  are  presented.  The  results  are  restricted  to  de¬ 
formation  for  which  isothermal  behavior  prevails  and  displacement  gra¬ 
dients  are  small.  Seven  types  of  mathematical  material  models  are  dis¬ 
cussed  that  can  be  used  to  describe  the  stress -strain -time  behavior  of 
physically  nonlinear  materials.  The  models  are  given  in  three- 
dimensional  tensorial  representation. 

The  energy -dissipation  properties  of  engineering  materials  are 
discussed  within  the  framework  of  the  above-mentioned  mathematical 
material  models,  and  the  damping  models  used  in  engineering  analyses  of 
vibratory  systems  are  summarized. 


Preceding  page  blank 
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DAMPING  CAPACITY  OF  SOIL  DURING  DYNAMIC  LOADING 


REVIEW  OF  MATHEMATICAL  MATERIAL  MODELS 

PART  I:  INTRODUCTION 

Background 

1.  A  great  number  of  engineering  materials  exhibit  hysteretie 
effects  under  cyclic  or  near-cyclic  loading  conditions,  even  when 
undergoing  small  deformation.  A-  prime  example  is  soil.  The  energy  - 
dissipation  characteristics  of  soil  can  be  observed,  for  example,  from 
the  attenuation  of  the  stress  pulse  during  one -dimensional  -wave  propa¬ 
gation  testa  or  from  the  formation  of  a  hysteretie  loop  during  a  load- 
unload  cycle  of  deformation.  The  hysteretie  character  of  the  stress - 
strain  behavior  of  soil  plays  an  important  role  in  the  determination  of 
the  response  of  soil  to  dynamic  loads  and  must  be  incorporated  in  the 
formulation  and  solution  of  soil  dynamics  boundary- value  problems. 

2.  The  term  damping  is  often  used,  in  a  general  sense,  to  de¬ 
scribe  the  energy -dissipation  properties  of  engineering  materials;  and 
various  mathematical  material  models  have  been  devised  to  represent 
such  phenomena.  The  approach  used  to  formulate  most  of  these  models  is 
based  on  the  methods  of  continuum  mechanics,  which  disregard  the  atom¬ 
istic  details  of  the  material  and  consider  its  gross  behavior  only. 

The  models  are,  therefore,  phenomenological  in  nature.  In  fact,  no 
satisfactory  explanation  of  damping  mechanism  on  the  atomistic  level  is 
available  at  the  present  time,  although  several  hypotheses  have  been 
set  f  >rth  in  recent  years.  For  the  solution  of  many  engineering  prob¬ 
lems,  however,  a  matneaatical  model  of  damping  that  is  based  on  the 
concepts  of  continuum  mechanics  is  quite  sufficient  and  useful.  Such  a 
model  must,  of  course,,  be  general  enough  to  describe  the  energy- 
dissipation  properties  of  the  material  under  various  spates  of  stress 
and  deformation. 

3.  To  date,  two  alternate  procedure  nave  been  employed  to 
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describe  the  hysteretic  behavior  or  the  energy -dissipation  properties 
of  soil.  Both  procedures  sre  based  on  the  methods  of  continuum  mechan¬ 
ics.  In  the  first  procedure,  jhe  hysteretic  behavior  of  soil  is  ac¬ 
counted  for  by  using  three  st.ts  of  time -independent  stress -strain  rela¬ 
tions:  one  set  for  initial  loading,  one  lor  unloading,  and  one  for 
subsequent  reloading.  The  mathematical  expressions  for  the  stress- 
strain  relations  are  usually  derived  from  cyclic  or  near-cyclic  experi¬ 
mental  data  by  curve -fitting  techniques.  A  set  of  criteria  or  logics 
is  also  specified  to  determine  whether  the  material  under  consideration 
is  being  loaded,  unloaded,  or  r  loaded  so  that  the  proper  stress -strain 
relations  can  be  used.  The  theory  of  plasticity  is  the  outcome  of  this 
procedure.  In  the  second  procedure,  a  time -dependent  mathematical 
model  is  postulated,  and  the  hysteretic  behavior  of  soil  is  accounted 
for  by  the  appearance  of  the  time  derivatives  (deformation  rate)  in  the 
model.  The  theory  of  viscoelasticity  is  the  outcome  of  this  procedure. 

4.  Both  of  these  techniques  have  been  used  extensively  in  the 
formulation  and  solution  of  soil  dynamics  boundary- value  problems.  The 
mechanism  of  energy  dissipation  in  the  first  procedure,  however,  is 
different  from  that  in  the  second  procedure.  At  the  present  time,  it 
is  not  clear  which  procedure  is  the  more  correct  one,  and  the  available 
experimental  information  on  the  stress -strain  behavior  of  soil  under 
dynamic  cyclic  loading  is  not  conclusive  enough  to  favor  one  or  the 
other.  For  reasons  of  practicability  or  mathematical  simplicity,  and 
depending  on  the  type  of  problem  being  considered,  some  investigators 
use  the  first  procedure  while  others  use  the  second.  A  subsequent  re¬ 
port  (Report  2)  will  review  current  laboratory  methods  of  determining 
the  damping  capacity  of  soil;  it  will  describe  specific  techniques  used 
to  evaluate  hysteretic  effects  within  the  framework  of  both  the  time- 
independent  and  the  time -dependent  mathematical  material  models. 

Purpose 

5.  The  purpose  of  this  report  is  to  review  and  document  various 
mathematical  material  models  that  can  be  used  to  describe  the 
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stress-strain-time  behavior  of  physically  nonlinear,  isotropic,  and 
homogeneous  materials.  Emphasis  is  placed  on  the  material  models  that 
exhibit  hysteretic  effects  under  cyclic  loading.  These  models  will 
serve  as  bases  for  tuture  work  for  the  evaluation  of  damping  capacity 
of  soils  under  dynamic  loads. 


Scope 

C.  Since  the  mathematical  material  models  presented  in  this  re¬ 
port  are  based  on  the  methods  of  continuum  mechanics,  a  brief  discus¬ 
sion  of  the  fundamental  concepts  of  the  theory  of  continuous  media  is 
given  in  Part  II.  Part  III  contains  a  summary  of  the  available  mathe¬ 
matical  material  models  in  their  general  functional  representation.  In 
Part  IV,  various  types  of  damping  models  commonly  used  in  engineering 
analyses  of  vibratory  systems  are  presented. 
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PAST  IX:  CONCEPTS  FROM  CONTINUUM  MECHANICS 

Basic  Principles 


7.  In  engineering  practice  it  is  convenient,  and  often  reason¬ 
able^  to  disregiv  ,  the  atomistic  details  of  »atter  and  consider  its 
gross  behavior  •  ,«ly.  It  is  also  customary  to,  adopt  the  hypothesis  that 
matter  can  be  replaced  by  a  mathematical  model  whose  kinematic  or  dy¬ 
namic  variables  are  piecewise  continuous  functions  of  spatial  coordi¬ 
nates  and  time,  Si’ch  a  medium  is  often  referred  to  as  a  continuum. 

The  motion  of  any  continuum  in  a  Galilean  frame  of  reference  is  gov¬ 
erned  by  the  following  laws: 

a.  Conservation  of  mass 

b.  Conservation  of  energy 

c.  Balance  of  linear  momentum 

d.  Balance  of  angular  momentum 

e.  Principle  of  inadmissibility  of  decreasing  entropy 
These  laws  constitute  the  basic  axioms  of  continuum  mechanics.1  If  me¬ 
chanical  energy  is  the  only  form  of  energy  to  be  considered  in  a  prob¬ 
lem,  the  above  principles  will  lead  to  the  continuity  equation* 


and  the  equations  of  motion 


0.,  .  +  p(F. 


a. )  -  0 
i ' 


(1) 

(2) 


where 

p  =  mass  density 

v^  =  components  of  velocity  vector 


*  Indices  take  on  values  1,  2,  or  3 •  A  repeated  index  is  to  be  summed 
out  over  its  range.  Comma  in  the  subscripts  represents  a  derivative. 
Quantities  are  referred  to  rectangular  Cartesian  coordinates  X.  . 
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o...  =  symmetrical  stress  tensor 
3-vl 

F^  =  components  of  body  force 

a..  =  components  of  acceleration  vector 

8.  Equations  1  and  2  constitute  four  equations  that  involve  ten 
unknown  functions  of  time  and  space:  the  mass  densit-"-  p  ,  the  three 
velocity  components  v^  ,  and  the  six  independent  stress  components 

ck  ..  .  The  body  force  components  F^  are  known  quantities  and  the  ac¬ 
celeration  components  ai  are  expressible  in  terms  of  the  velocity 
components  vi  .  Obviously  equations  1  and  2  are  inadequate  to  deter¬ 
mine  the  motion  or  deformation  of  a  medium  when  subjected  to  external 
disturbances  such  as  surface  forces  and/or  displacements.  Therefore, 
six  additional  equations  relating  the  ten  unknown  variables  p  ,  v^  , 
and  are  required.  L:  the  field  of  continuum  mechanics,  such  re¬ 

lations  art-  stated  by  constitutive  equations,  which  relate  the  stress 
tensor  o..  uo  the  deformation  and  time  rate  of  deformation.  The  dif- 
ference  between  the  constitutive  equations  and  the  field  equations 
(equations  1  and  2)  is  that  the  latter  are  applicable  to  all  materials, 
whereas  the  former  represent  the  intrinsic  response  of  a  particular  ma¬ 
terial.  Furthermore,  a  constitutive  equation  provides  a  mathematical 
description  or  definition  of  an  ideal  material  rather  than  a  statement 
of  a  universal  law. 

9.  The  general  form  of  a  constitutive  equation  for  isothermal 
conditions  may  be  expressed  by  the  functional  form 


’  0 


(3) 


relating  the  ten  unknown  variables  p  ,  v^  ,  and  .  Equation  3 
can  be  written  in  a  more  definite  form  relating  the  stress  tensor  to 
the  strain  tensor  and  the  deformation-rate  tensor,  i.e.. 


a.  .  =  f . ,(e  ,d  ) 
ij  lj'  mn’  rs/ 


(4) 


where 


e  =  strain  tensor 
mn 

d  =  deformation -rate  tensor 
rs 
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If  the  displacement  gradients  are  small,  the  strain  tensor  can  be  ex¬ 
pressed  in  terms  of  the  components  of  the  displacement  vector  vu  by 
the  following  relation 


a 

~mn 


-  h  (u 
2  v  m,n 


+  u  ) 


(5) 


The  deformation -rate  tensor  is  related  to  the  components  of  the  veloc¬ 
ity  vector  vi  ,  i.e. 


drs  =  I  ^vr,s 


+  V  ) 

e>  v»  • 


s,r' 


(6; 


The  mass  density  p  can  be  related  to  dys  or  through  equations 

1,  5,  and  6.  Equations  1,  2,  and  4,  therefore,  constitute  ten  equa¬ 
tions  in  ten  unknowns  and  will  lead,  ifr  conjunction  with  the  kinematic 
relations  given  by  equations  5  and  6,  to  a  complete  description  of  the 
boundary -value  problem.  In  addition  to  the  above-mentioned  equations, 
boundary  conditions  in  terms  of  boundary  displacements  and/or  surface 
tractions  must  also  be  specified  to  define  the  particular  problem  of 
interest. 


Requirements  for  Constitutive  Equations 


10.  In  order  for  the  constitutive  equation  (equation  4)  to  de¬ 
scribe  a  physical  material  adequately,  the  response  function  f. .  must 
be  form  invariant  with  respect  to  rigid  motion  of  spatial  coordinates. 
This  requirement  terns  from  the  fact  that  the  response  of  a  material  is 
independent  of  the  motion  of  the  observer.  Furthermore,  must  be 
expressed  in  tensor  language  to  ensure  that  the  constitutive  equation 
is  invariant  to  coordinate  transformations.  In  addition  to  the  above- 
mentioned  invariance  principles,  the  response  function  f^..  must  be 
consistent  with  the  general  principles  of  conservation  or  balance  of 
mass,  momentum,  and  energy.  If  it  is  assumed  that  the  material  under 
consideration  is  isotropic  and  elastic  (linear  or  nonlinear),  every 
principal  axis  of  strain  must  also  be  a  principal  axis  of  stress,  i.e., 


t 


■% 


if 


€12  "  e23  "  631  "  °  * 


ID 


must  show  that 


°12  "  °23  °31  ~  ° 


Also,  to  avoid,  stress -induced  anisotropy  for  an  isotropic  elastic 
medium,  f^  must  ensure  that  a  zero  state  of  strain  corresponds  to  a 
zero  or  a  specified  scalar  state  of  stress. 


The  General  Form  of  Isotropic  Constitutive  Equation 


11.  The  general  form  of  the  response  function  f . ,  satisfying 

the  invariance  principles  and  the  requirements  of  isotropy  was  derived 

2 

by  Rivlin  and  Ericksen.  Functions  satisfying  these  conditions  are 
called  hemitropic  functions  of  their  arguments.  A  hemitropic  poly¬ 
nomial  F  of  two  symmetric  second-rank  tensor  variables  A  and  B 
admits  a  representation  of  the  form 


F(A,B)  =  fiQI  +  fi. jA  +  A2  +  B  +  ^B2  +  fi^AB  +  BA.)  +  ^(A2B  +  BA2) 

+  ^?(AB2  +  b2a)  +  ^8(a2b2  +  b2a2)  (7) 

in  which  I  =  6^  =  the  Kronecker  delta  ;  and  the  coefficients  0^  ,  . . .  , 

fin  are  scalar-valued  functions  of  the  ten  joint  invariants  of  A  and 
0  3 

B  -  The  ten  joint  invariants  of  A  and  B  are  given  by 


tr(A) 

tr(A2) 

tr(A3) 

tr(B) 

tr(B2) 

tr(B3) 

tr(AB) 

tr(A2B) 

tr(B2A) 

tr(A2B2) 

where  tr  =  trace  of  ,  indicating  the  sum  of  the  diagonal  terms  of  a 
square  matrix. 

12.  According  to  equation  7>  the  admissible  form  of  the  isother¬ 
mal  constitutive  equation  (equation  4),  or  f . ,  ,  takes  the  following 

1 J 

form 
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cr.  =  tia5.  .  +  tJt  e .  .  +  n0e.  e  ;  +  q~d.  .  +  r>,d.  d  .  +  n,-(e.  d  .  +  d.  e  .) 
10  '0  10  '1  10  '2  Ml  mo  '2  13  4  im  mj  '5  im  mo  im  mj' 

+  tjz-( e .  e  d  .  +  d.  e  e  .)  +  n_(e.  d  d  .  +  d.  d  e  ) 

6V  im  mn  no  im  mn  no  ]7K  im  mn  no  im  mn  nj' 

+  n0(£-  e  d  d  .  +  d.  d  e  e  . )  fc 

‘8  s  im  mn  np  po  im  mu  np  pj 


The  response  coefficients  ,  ...  ,  rig  are  scalar-valued  functions  of 
the  ten  joint  invariants  of  the  strain  tensor  and  the  deformation -rate 
tensor.  In  view  of  equations  5,  6,  and  8,  the  ten  joint  invariants  of 
the  strain  tensor  and  the  deformation -rate  tensor  take  the  following 
representation: 


I.  = 

€ 

• 

I. 

=  d 

1 

SS 

1 

ss 

I2  = 

e.  e  , 
ts  st 

• 

I 

2 

^ts^st 

T  = 

a  e  , 

• 

I„ 

=  d  d  d 

3 

ts  sr  rt 

3 

ts  sr  rt 

M  = 

6ts^st 

R  = 

Ss^sr^rt 

N  = 

etsesrdrt 

Q  = 

6tseskdkr<i] 

These  response  coefficients  t»-ke  various  forms  for  different  materials 
and  must  he  determined  from  experimental  observation.  There  is,  how¬ 
ever,  no  a  priori  reason  for  requiring  that  all  response  coefficients 
appear  in  the  constitutive  equations  for  all  materials.  Some  of  the 
response  coefficients  may  vanish  for  some  materials.  Accordingly,  a 
constitutive  equation  with  two  response  coefficients,  such  as  Hooke's 
law,  is  as  valid  and  significant  as  a  constitutive  equation  that  in¬ 
cludes  the  nine  possible  response  coefficients  t|q  tjg  ,  The  dif¬ 

ference  between  the  two  equations  is  in  their  range  of  application. 

13.  It  should  be  pointed  out  that  equation  9  is  not  based  on 
thermodynamic  considerations  and  the  response  coefficients  are  not,  in 
general,  related  to  a  single  potential  function.  A  special  form  of 
equation  9  based  on  thermodynamic  considerations  was  also  developed 
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by  Schapery  and  has  been  used  extensively  for  stress  analysis  of  visco^ 
elastic  materials . 

14.  Presently  there  are  several  forms  of  isotr  1c  constitutive 
equations  available  that  can  be  used  to  describe  the  stress- strain- time 
response  of  various  materials.  They  are  all  derived  from  equation  7 
or  9  and  are  presented  in  the  following  part. 
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PART  III:  REVIEW  OF  ISOTROPIC  CONSTITUTIVE  EQUATIONS 

'  General 


15.  The  theoretical  foundation  of  isotropic  constitutive  equa¬ 
tions  is  based  on  the  Rivlin-Ericksen  equation  presented  in  Part  II  (equa^ 
tion  7)-  Various  classes  of  constitutive  equations  have  been  developed 
in  recent  years  utilizing  equation  7  or  9*  A  summary  of  these  equations 
is  presented  in  the  following  sections. 


Elastic  Materials 

16.  For  elastic  materials,  the  state  of  stress  is  a  function  of 
the  current  state  of  strain,  i.e. 


aij  =  fi^6mn^ 


The  response  function  f^  has  the  same  form  as  equation  7  where 
A  =  emn  and  B  =  0  .  Accordingly,  the  constitutive  equation  for  elas¬ 
tic  materials  becomes 


Equation  12  is  often  referred  to  as  the  Cauchy  elastic  constitutive 
equation.1  The  response  coefficients  0O  ,  0^  ,  and  0g  are  func¬ 

tions  of  the  three  strain  invariants  1^  ,  Ig  ,  and  I^  given  by  equa¬ 
tion  10.  Equation  12  can  describe  the  mechanical  behavior  of  various 
types  of  elastic  materials  by  properly  selecting  the  forms  of  0^  , 

01  ,  and  0g  .  For  example,  if  0g  =  0  ,  0^  =  2 p  ,  and  0Q  =  XI  , 
equation  12  degenerates  to  Hooke's  law,  where  n  and  X  are  the  Lame 
constants.  For  a  second -order  elastic  stress -strain  relation,  the  re¬ 
sponse  coefficients  take  the  forms  0O  =  D  ,  0.  =  2\x  +  CL. ,  and 
0q  =  XIX  +  AI£  +  BIg  where  A  ,  B  ,  C  ,  and  D  are  material  con¬ 
stants.  Therefore,  a  second-order  elastic  stress -strain  law,  based  on 
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equation  12,  involves  six  material  constants.  Paralleling  the  above 
procedure,  higher  order  and  more  complicated  elastic  stress -strain  re¬ 
lations  can  be  developed. 

17.  The  counterpart  of  equation  12,  which  is  derived  from  the 
conservation  of  energy  and  is  referred  to  as  the  hyperelastic  constitu¬ 
tive  equation,  is  given  as 


(13) 


5 

where  U  =  U(l^,Ig,l^)  is  the  strain-energy  density  function.  Com¬ 
parison  of  equation  13  with  equation  12  indicates  that  hyperelastic 
materials  are  special  forms  of  the  Cauchy  elastic  materials,  where 


d  -  5u 


O'  = 


dU  .  d  _ 


,  dU 


(14) 


Consequently,  a  Cauchy  elastic  material  is  hyperelastic  if  the  response 
coefficients  0^  ,  0^  ,  and  02  are  related  in  the  following  manner 


For  linear  elastic  materials 


*01 

(15a) 

^1  2 

(15b) 

*2  »0 

(15c) 

■  H  *  ^ 

(16) 

ar.d  both  equations  12  and  13  yield  the  same  results.  The  effect  of  the 
thermodynamic  restrictions,  therefore,  is  not  evident  when  considering 
linear  elastic  materials.  However,  the  effect  of  thermodynamic 
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restrictions  becomes  very  pronounced  when  considering  higher  order 
stress -strain  laws.  For  example,  a  second-order  elastic  stress -strain 
law  formulated  from  the  Cauchy  elastic  constitutive  equation  involves, 
as  was  shown  previously,  six  material  constants.  On  the  other  hand,  a 
second-order  elastic  stress -strain  law  based  on  the  hyperelastic  con¬ 
stitutive  equation  (equation  13 )  involves  only  five  material  con- 
stants,  since  the  strain  energy  function  for  second -or^er  hy\.erelastic 
materials  takes  the  form  U  =  g-  1^  +  +  C.^  where 

,  and  are  material  constants.  The  effect  of  thermodynamic 

restrictions,  in  this  case,  is  to  reduce  the  number  of  material  con¬ 
stants  from  six  to  five.  However,  the  physical  implications  of  this 
reduction  in  material  constants  is  not  clear  at  the  present  time. 

18.  The  inverse  forms  of  the  Cauchy  elastic  and  hyperelastic  con¬ 
stitutive  equations,  resulting  in  strain-stress  laws,  are  given  as 


e«  ‘  eoftia 


+  Via 


0„<J.  0  , 

2  im  mj 


(17) 


for  the  Cauchy  elastic  materials,  and 


e.  . 

U 


_  cty 


6.  . 

id 


+  2 


5^  °i3  +  3  Sjj  W) 


(18) 


for  hyperelastic  materials.  The  response  coefficients  0Q  ,  ,  and 

02  ,  as  well  as  the  complementary-energy  density  function  ^  ,  are 
functions  of  the  following  stress  invariants 


J1  =  ass 

J2  =  atsast  (19) 

J3  ~  atsCTsrart 


For  linear  elastic  materials,  the  complementary-energy  density  function 
is  given  as 


(20) 


t 


v 


♦  = 


1  +  V 

2E 


v  t2 
2E  °1 


where  E  and  v  are  Young's  modulus  of  elasticity -and  Poisson's  ratio, 
respectively;  If  equation  17  is  utilized  to  derive  the  strain-stress 
law  for  linear  elastic  materials,  the  response  coefficients  take  the 
forms  02  =  0  ,  0.^  =  — ,,  and  0Q  =  “  e  J1  '  ASain>  as  expected, 

both  the  Cauchy  elastic  and  hyperelastic  constitutive  equations  yield 
the  same  strain-stress  laws  for  linear  elastic  materials. 


Hypoelastic  Materials 

7 

19*  The  theory  of  hypoelasticity  was  formulated  by  Truesdell. 

The  theory  is  formulated  in  terms  of  rates  of  both  stress  and  deforma¬ 
tion  and  is  intended  to  describe  the  mechanical  behavior  of  path- 
dependent  materials.  As  was  pointed  out  in  the  previous  section,  the 
state  of  stress  for  an  elastic  material  is  a  function  of  the  current 
state  of  strain  and  is  independent  of  the  path  followed  to  reach  that 
state.  For  real  materials,  earth  materials  in  particular,  the  final 
state  of  stress  is  a  function  of  the  final  state  of  strain  as  well  as 
the  stress  path  used  to  reach  the  final  state.  Hypoelastic  theory  pre¬ 
dicts  this  type  of  behavior  and  may  be  used  to  mode,!  the  stress -strain 
behavior  of  soil. 

20.  The  basic  constitutive  equation  of  hypoelastic  materials  is 
expressed  in  the  form 


s  =  f(d,s) 


(21) 


.where 


A 
0 . 


s  =  nondimensional  stress  flux  tensor  = 

2  \i 

.Si 

2n 


s  =  nondimensional  stress  tensor 
0^  =  stress  flux  tensor 


d  =  deformation -rate  tensor  (equation  6) 

1 

Jaumann’s  form  of  stress  flux  tensor  is  given  by 
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(22) 


\ 


bo.  . 


"  "5t  *  °«,kTk  +  °ik“k3  -  “kAk 


where 


id  =  o  (v  -  v 
rs  2  v  r,s  s 


r) 


(93) 


is  the  Eulerian  spin  tensor.  This  definition  of  stress  flux  is  not 

unique,  and  various  other  forms  developed  hy  Truesdell  and  Oldroyd  are 

i 

available,  which  differ  from  Jaumann’s  form  in  the  terms  containing  the 
spin  tensor.  However,  at  the  present  time  there  is  no  a  priori  reason 
for  utilizing  any  specific  form  of  the  stress  flux  tensor;  additional 
physical  postulates,  or  experimental  information,  are  required  in  order 
to  determine  a  preferred  form  of  this  tensor.  For  static  problems  or 
for  dynamic  -problems  with  nearly  irrobational  displacement  fields,  the 
stress -rate  tensor  (the  first  term  on  the  right  side  of  equation  22) 
may  be  used  instead  of  the  stress  flux  tensor. 

21.  The  hypoelastic  response  function  f(d,s)  is  of  the  same 
form  as  equation  7  where  A  =  d  and  B  =  s  .  Accordingly,  the  hypo- 
elastic  constitutive  equation  (equation  21 )  becomes 

2  2 

s  =  aQI  +  ctjd.  +  a2d  +  o^s  +  a^s  +  (ds  +  sd) 

+  ag(d2s  +  sd^)  +  (ds2  +  s^d)  +  (d2s2  +  scd^)  (24) 

where  cCq  ,  ...  ,  Gg  are  scalar-valued  functions  of  the  ten  joint  in¬ 
variants  of  d  and  s  (joint  invariants  of  d  and  s  are  obtained 
from  equation  8  by  substituting  d  and  s  for  A  and  B  ,  respec¬ 
tively).  The  response  coefficients  oj-,  ,  ...  ,  are  further  re¬ 
stricted  by  the  first  hypothesis  of  hypoelasticity:  "No  constitutive 
coefficients  of  a  hypoelastic  material  shall  carry  a  dimension  inde¬ 
pendent  of  the  dimension  of  stress."  The  consequence  of  this  hypothesis 
on  the  hypoelastic  constitutive  equation  is  that  all  terms  containing 
second  and  higher  powers  of  d  must  vanish.  Thus,  the  response  coef¬ 
ficients  ag  =  =  ag  =  0  5  ,  and  must  be  independent  of 

d  anci.  functions  of  s  alone;  and  o^q  >  ,  and  must  be  of  degree 

one  in  d  .  Imposing  the  above  restrictions  on  the  response  coefficients 
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in  equation  2k,  the  hypoelastic  constitutive  equation  reduces  to 


p 

§  =  ijP0i  +  pxa  +  i1P2§  +  ep3i  +  -f  (ds  +  sd)  +  s2 


p, 


+  EpgS  +  Fp?I  +  (ds2  +  S2d)  +  EP^*  +  FP10s  +  Fp-^s*  (25) 


-2 


-2 


in  which  Pq  ,  ...  ,  P  are  dimensionless  functions  of  the  three  prin¬ 
cipal  invariants  of  s  only,  and  E  and  F  are  defined  as 


(26) 

F  =  s  s  d 
mn  np  pm 

22.  From  equation  25,  it  is  apparent  that  the  constitutive  equa¬ 
tions  of  hypoelasticity  are  coupled  differential  equations  of  first  or¬ 
der.  Furthermore,  the  differential  equations  are  homogeneous  in  time. 
To  obtain  a  unique  solution  to  these  equations,  some  initial  conditions 
that  are  consistent  with  the  invariance  principles  must  be  prescribed. 
The  integration  of  the  differential  equations,  for  a  given  stress  path, 
leads  to  stress-strain  relations.  Thus,,  a  relation  between  stress  and 

strain  is  the  outcome  of  the  theory. 

7 

23.  Truesdell  has  defined  various  classes  of  hypoelastic  mate- 

» 

rials  that  are  characterized  by  the  highest  degree  of  s  appearing  in 
equation  25.  If  the  right  side  of  equation  25  is  independent  of  ’s  , 
the  material  is  called  hypoelastic  material  of  grade  zero.  In  this 
case  equation  25  reduces  to 


/\ 

s 


Vo1 


+  P-jd 


(27) 


which  .is  directly  similar  to  the  constitutive  equation  for  a  linear, 
isotropic,  elastic  material  if  PQ  =  Jjj  and  p^  =  1  ,  i.e. 


10 


(28) 
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( 


If  the  right  side  of  equation  25  contains  up  to  the  first  power  of 
stress,  the  material  is  called  hypoelastic  of  grade  one  and  so  on. 

24.  In  establishing  the  constitutive  equation  of  hypoelastic 
materials,  no  assumptions  were  made  in  regard  to  the  magnitude  of 
stresses  and  strains  involved;  hence,  the  theory  is  applicable  for  all 
motions . 


Hygrosteric  Materials 

25.  The  term  hygrosteric  is  used  for  materials  having  a  consti¬ 
tutive  equation  of  the  form 


s  =  f(d,s,p)  (29) 

where  the  response  function  f  is  a  polynomial  In  d  and  s  with  co¬ 
efficients  depending  on  the  density  p  .  Due  to  the  scalar  character 
of  the  density  p  ,  the  formulation  of  constit;  rve  equations  for  hygro¬ 
steric  materials  is  the  same  as  for  the  hypoelastic  materials  except 
that  the  response  coefficients  appearing  'in  equation  25  may  be  con¬ 
sidered  functions  of  p  or  the  dimensionless  ratio  p/pQ  where  pQ 
is  the  initial  mass  density  of  the  raterial. 

Stokesian  Fluids 


26.  Stokesian  fluids  are  characterized  by  the  constitutive  equa¬ 
tions  of  the  following  form 


a..  =  f..(d  ) 
ij'  rs 


(30) 


The  state  of  stress,  therefore,  is  a  function  of  the  current  rate  of 
deformation.  In  view  of  equation  9 5  equation  30  becomes 


(31) 
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■where  q 


0 


n3  >  and 


are  functions  of  the  three  deformation-rate 


invariants  given  by  equation  10.  Equation  31  can  describe  the  behavior 
of  various  fluids  by  proper  selection  of  the  response-  coefficients,  n  , 


,  and 
where  jo. 


=  -P  +  X 


\  •  example,  if  tIu  »  0  ^  =  2^  ,  and  r\Q  -  ±  . 

v  and  X^  are  the  shear  and  dilatati  «,1  viscosity  coefficients, 
respectively,  and  P  is  pressure,  equation  31  degenerates  to  the  con¬ 
stitutive  equation  of  linear  viscous  fluids,',  i.e. 


B«  ’  -rei,1  +  Vl8i:,  +  Via 


10 


iO 


(Si.) 


Viscoelastic  Materials 


27.  The  constitutive  equation  of  viscoelastic  materials  is  ex¬ 
pressed  by 


(e  ,d  ) 
'  mn’  rs' 


(33) 


where  the  isotropic  function  fy  is  given  by  equation  9.  Equation  33 
reduces  to  the  Cauchy  elastic  constitutive  equation  (equation  12)  if  de¬ 
pendence  on  dys  disappears.  Various  classes  of  viscoelastic  mate-* 
rials  can  be  described  by  equation  33  (or  equation  9)  by  proper  se?  c- 
tion  of  the  response  coefficients  Pq  ,  ...  ,  Pg  .  For  example,  in' 

P0  ~  +  Vl  ’  ^1  =  ^  }  ^3  =  2^v  *  an<^  ^2  ~  ^4  ~  ^5  ~  ^6  ~  ^r( 

=  i\q  -  0  ,  equation  9  reduces  to  the  constitutive  equation  of  the 
Kelvin -Vpigt  material,  i.e. 


=  XIlsiJ  *  WiJ  +  Vd  +  Vid  (34) 

Equation  34  is  the  three-dimensional  representation  of  the  parallel 

spring-dashpot  model  used  in  the  theory  of  vibration.  Taxing  tj  =  V 

5  o 

=  Py  =  hg  -  0  in  equation  9,  a  second-order  viscoelastic  constitutive 
model,  often  referred  to  as  the  nonlinear  Kelvin  solid,  can  result,  i.e. 
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•H 


(35) 


°id 


Oi; 


1  id 


c.  iiii  uy 


Equation  35  includes  tensorial  as  well  as  scalar  nonlinearities  and  can 
no  longer  be  presented  by  a  simple  parallel  spring  and  dashpot  model. 
Many  complex  aspects  of  the  physical  behavior  of  real  materials,  such 
as  shear -dilatancy  phenomenon  and  the  Poynting  effect,'1'  can  be  de¬ 
scribed  by  equation  35- 

28.  Viscoelastic  constitutive  equations  can  also  be  formulated 
in  integral  forms,  e.g. 


o. . 

id 


t 

■/ 


dl, 

X  ( 11  "  T  )  dT  8ijdT  +  2  /  ^  ( t  "  T ) 


de. 

dx 


dx 


(36) 


where  X  and  p  are  referred  to  as  relaxation  moduli  or  memory  func¬ 
tions  and  (t  -  x)  is  an  elapsed  time.  Other  forms  of  viscoelastic  con¬ 
stitutive  equations,  e.g.,  series  forms  with  differential  operators  as 
arguments,  are  also  available  that  can  be  used  to  characterize  various 
viscoelastic  materials.1  Viscoelastic  constitutive  equations  that  are 
expressed  in  the  form  of  equation  9  (e.g.,  equations  34  and  35) ,  how¬ 
ever,  possess  a  certain  mathematical  simplicity  that  makes  them  more 
attractive  for  engineering  analyses  than  other  forms  of  viscoelastic 
constitutive  equations. 


Viscoplastic  Materials 

29.  To  describe  the  mechanical  behavior  of  rate -dependent  mate¬ 
rials  that  are  corapaetible,  i.e.  exhibit  time -independent  as  well  as 
time -dependent  hysteretic  effects,  constitutive  equations  of  the  fol¬ 
lowing  form  are  often  used: 


a. .  =  f. .  (a  ,e  ,d  ) 
id  '  pq  mn’  rs' 


(37) 


where  0  are  the  components  of  the  stress -rate  tensor, 
tropic  response  function  f . .  is  a  hemitropic  function  of 

10  Q 

and  d  and  takes  the  following  representation: 


The  iso- 
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a.  .  =  A.6. . 
10  0  13 


a,  a.  . 
1  10 


2  im  m3 


+  A„(a.  e  .  +  €.  cr  .)  +  Ao (d.  d  .  +  d.  a  .)  +  AAd.  e  a  . 
7  im  mj  xm  m3  o  im  mo  im  mo  9  im  mn  no 


+  A,a£.  d  s  .  +  A-.d.  d  d  .  +  A.Ad.  &  d  . 
10  im  mn  no  11  im  mn  no  12  xm  mn  no 


m 


The  response  coefficients  Aq  ,  ...  ,  are  scalar  functions  of  fif¬ 

teen  3‘oint  invariants  of  cr  ,  e  „  ,  and  d  .  Equation  38  reduces 

pq  mn  rs 

to  the  Cauchy  elastic  constitutive  equation  (equation  12)  if  dependence 
on  d  and  d  disappears.  If  dependence  on  d  and  e  dis- 
appears,  equation  38  will  reduce  to  the  constitutive  equation  of  the 
Stokesian  fluid  (equation  31)*  If  dependence  on  d  disappears,  equa¬ 
tion  38  will  reduce  to  the  constitutive  equation  of  the  nonlinear  Kel¬ 
vin  solid  (equation  35).  Equation  38  is,  therefore,  very  general,  and 
by  proper  selection  of  the  response  coefficients  AQ  ,  ...  ,  A^g  ,  it 
can  be  used  to  describe  the  mechanical  behavior  of  many  complex 
materials . 


Summary  of  Constitutive  Equations 

30’.  In  summary,  seven  basic  types  of  isotropic  constitutive  equa 
tions  that  can  be  utilized  to  describe  the  mechanical  behavior  of  real 
materials  have  been  discussed.  These  constitutive  equations  are  all 
expressed  in  general  functional  forms.  The  specific  forms  of  the  func¬ 
tionals,  for  any  particular  material,  must  be  determined  through  exper¬ 
imental  observation  or  by  physical  postulates.  The  constitutive  equa¬ 
tions  are: 

a.  Cauchy  elastic  materials  (equation  12) 

b.  Hyperelastic  materials  (equation  13) 

c.  Hypoelastic  materials  (equation  25) 

d.  Hygrosteric  materials  (equation  29) 

e.  Stokesian  fluids  (equation  31) 

f.  Viscoelastic  materials  (equation  9) 

g.  Viscoplastic  materials  (equation  38) 
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Equations  12  and  13  are  used  to  describe  the  behavior  of  elastic  mate¬ 
rials  (linear  or  nonlinear).  Equations  25  and  29  are  used  to  describe 
the  behavior  of  rate -independent  hysteretic  (or  plastic)  materials. 
Equation  31  is  used  to  describe  the  behavior  of  various  fluids.  Equa¬ 
tion  9  is  used  to  describe  the  behavior  of  various  fate -dependent  (visco¬ 
elastic)  materials.  Equation  38  is  used  to  describe  the  behavior  of 
rate -dependent  hysteretic  (viscoplastic)  materials. 

21.  In  the  case  of  hysteretic  materials ,  the  constitutive  equa¬ 
tions  insist  be  accompanied  by  a  set  of  criteria  (or  logics)  defining 
initial  loading,  unloading,  and  subsequent  reloading  of  the  material. 
These  criteria  must  be  stated  in  terms  of  the  invariants  of  the  stress 
and/or  strain  tensors  in  order  to  remain  invariant  with  respect  to 
rigid  motion  of  spatial  coordinates.  The  hysteretic  behavior  of  the 
material  is  then  taken  into  account  by  using  different  sets  of  consti¬ 
tutive  equations  for  loading,  unloading,  and  reloading.  There  are  sev¬ 
eral  criteria  for  defining  loading,  unloading,  and  reloading  that  are 
presently  being  used  in  the  solution  of  two-dimensional  boundary -value 
problems  for  hysteretic  materials.  These  criteria  are  not  necessarily 
conclusive  and  although  they  yield  satisfactory  results  for  certain 
well-deiined  stress  paths,  they  might  lead  to  unrealistic  material  be¬ 
havior  under  more  complex  loading  conditions.  A  summary  of  these  cri¬ 
teria  is  given  in  the  following  section. 

Criteria  for  Initial  Loading,  Unloading,  and  Reloading 

32.  One  common  criterion  for  differentiating  among  initial  load¬ 
ing,  unloading,  and  reloading  is  based  on  the  time-rate  (or  increment) 

of  the  first  invariant  of  the  stress  tensor  0  .  According  to  this 

ss 

criterion,  0  >  0  defines  loading  and  0  <0  defines  unloading, 

s  s  ss 

Whenever  0  >0  but  c  is  less  than  its  previous  maximum  value,  the 

ss  ss 

material  is  assumed  to  be  reloading.  A  similar  criterion  based  on  the 

time-rate  of  the  first  invariant  of  the  strain  tensor  e  can  also  be 

ss 

used  to  define  various  loading  conditions.  Both  of  these  criteria,  how¬ 
ever,  will  lead  to  unrealistic  material  behavior  during  a  load-unload 
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fT' 


cycle  where  the  medium  simultaneously  experiences  loading  in  shear  and 
unloading  in  pressure  or  vice  versa. 

33*  Another  criterion  is  based  on  the  time-rate  of  the  octa¬ 


hedral  shearing  stress  t  ^  .  According  to  this  criterion,  t  ^  >  0 
defines  loading  and  t  ^  <  0  defines  unloading.  The  material  is  as¬ 
sumed  to  be  reloading  whenever  tqc^  >  0  and  t  .j.  is  less  than  its 
previous  maximum  value.  The  condition  tqc^.  =  0  is  referred  to  as  a 


neutral  state  of  loading. 


The  t  ,  criterion  is  less  controversial 
oct 


than  the  6  or  e  criterion, 
ss  ss 


34*  In  an  attempt  to  overcome  some  of  the  theoretical  and  exper¬ 


imental  problems  that  can  arise  through  the  use  of  the  t  ^  or  cfg 
criterion,  a  combination  of  the  two  is  often  utilized  that  separates 


the  deformation  into  the  hydrostatic  and  the  deviatoric  components. 

The  or  criterion  is  used  for  the  hydrostatic  part  of  the  de forma- 

tion,  and  the  t  .j.  criterion  is  used  for  the  deviatoric  components. 

In  this  manner,  it  is  possible  for  the  material  to  unload  in  shear 

while  loading  in  pressure  or  vice  versa. 

35.  The  rate  of  work  a.  .k.  .  is  also  used  as  a  criterion  for 

ij  ij 

defining  various  loading  conditions  in  a  cyclic  test.  According  to 


this  criterion,  a.  .e,  .  >0  defines  loading  and  reloading;  a.  .e.  .  <0 

defines  unloading;  and  the  condition  a.  ,6.  .  =  0  is  referred  to  as  a 

**•  J  J 

neutral  state  of  loading.  'The  rate  of  work  criterion  and  the  tQc^ 
criterion  are  used  extensively  in  the  theory  of  plasticity,  and  they  are 


essentially  the  same  if  the  constitutive  equation  associated  with  them 


is  linear  from  the  tensorial  point  of  view. 

36.  The  neutral  states  of  loading  associated  with  the  toc^.  or 
a. .e. .  criterion  impose  certain  restrictions  on  the  material  constants 

10  10 

in  the  constitutive  equations  for  loading  and  unloading  and  require 
special  considerations.  The  material  constants  must  be  chosen  so  that 
the  loading  a.nd  unloading  constitutive  equations  become  identical  when¬ 
ever  tocj.  =  0  or  =  0  >  i.e.  neutral  loading.  This  require¬ 

ment  must  be  met  in  order  to  obtain  a  unique  solution  for  a  given 


boundary -value  problem  involving  cyclic  loading  conditions. 
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PART  IV:  SUMMARY  OF  DAMPING  MODELS  USED  IN 
ENGINEERING  ANALYSES  OF  VIBRATORY  SYSTEMS 


t 


37.  Most  of  the  damping  models  used  in  engineering  analyses  of 
vibratory  systems  stem  from  one -dimensional  representation  of  the  iso¬ 
tropic  constitutive  equations  presented  in  Part  III  and,  therefore,  are 
limited  in  their  application.  The  models  can  be  divided  into  two  gen¬ 
eral  categories:  viscous  damping  models  and  hysteretic  damping  models. 
Various  forms  of  these  models  are  discussed  in  the  following  sections 
and  are  given  in  their  three-dimensional  representation  for  a  broader 
application  and  utilization.  Appendix  A  gives  a  bibliography  of  source 
material  related  to  the  model  category  and  form. 


Viscous  Damping  Models 

38.  The  viscous  damping  models  are  frequently  represented  in 
terms  of  rheological  models  consisting  of  linear  springs  and  viscous 
elements  or  dashpots.  The  viscous  elements  account  for  the  energy- 
dissipation  properties  of  the  material,  while  the  springs  represent  the 
elastic  behavior  of  the  material.  Three  types  of  such  models  are  com¬ 
monly  used  in  engineering:  Kelvin-Voigt,  Maxwell,  and  standard-linear 
solid  models.  The  Kelvin-Voigt  model  has  be  in  used  more  extensively 
due  to  its  mathematical  simplicity.  The  three  models  are  presented  in 
the  order  of  increasing  complexity. 

Kelvin-Voigt  model 

39*  The  rheological  model  of  Kelvin-Voigt  h.atorial  consists  of  a 
linear  spring  and  a  dashpot  in  parallel  and  corresponds  to  a  one -degree - 
of -freedom  damped  system  in  the  theory  of  vibrations.  The  constitutive 
equation  of  Kelvin-Voigt  material  is  given  by 


a.  . 

ij 


=  (X^  + 


X  L) 

v  1 


6.  . 


2ue. .  +  2u  d. . 
10  *v  1J 


(39) 


Equation  39  involves  four  material  constants,  X  ,  X^  ,  p  ,  and  p^  , 
that  must  be  determined  experimentally.  The  one -dimensional 
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representation  of  equation  39?  which  is  most  commonly  used  in  engineer¬ 


ing  application,  can  be  obtained  by  allowing  =  €' 


23  =  S22 


=  =  ^22  =  d!3  =  ^3  =  ^22  =  ^33  =  ®  *  Accordingly,  in  view  of 

equations  5?  6,  and  10,  the  one -dimensional  representation  of  equation 
39  becomes 


a. ,  =  Xe.,  -  +  X 
11  11  v  dt 


+  2'i6n  +  2^v  "dF 


(40a) 


=  Xe..  +  X 

11  v  dt 


(40b) 


where 


is  the  strain  in  the  direction  of  motion  in  a  uniaxial 


strain  configuration  and  a ^  the  lateral  stress  required  to  pre¬ 
vent  lateral  strain.  Equation  40a  can  be  written  in  a  more  compact 
form  by  collecting  terms,  i.e. 


all  =  <X  +  2^6H  +  (\  +  V 


The  terms  X  +  2p  and  X^  +  2py  correspond  to  the  elastic  modulus  M 
and  the  distributed  viscosity  q  ,  respectively,  associated  with  the 
conditions  of  uniaxial  strain,  In  view  of  the  above  terminology,  equa¬ 
tion  41  becomes 


°u  ’  M£1X  +  1 


Equation  42  is  the  counterpart  of  the  differential  equation  of  motion 
for  the  free  vibration  of  a  one -degree -of -freedom  damped  system  with 
equivalent  spring  constant  and  damping  coefficient  cg^  given  by 


i  _  Ma 
Xeq  ”  h 


(43a) 


c  =  32. 
eq  h 


(43b) 


where  a  and  h  are  area  and  length,  respectively.  In  practice,  a 
corresponds  to  the  cross-sectional  area  of  a  continuous  system 
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(a  column)  and  h  corresponds  to  the  grid  spacing  of  a  discrete  system, 
(a  lumped -parameter  model)  that  is  mathematically  equivalent  to  the  con¬ 
tinuous  system. 

40.  Various  nonlinear  forms  of  Kelvin -Voigt  constitutive  equa¬ 
tions  are  also  available  that  can  be  used  to  describe  the  stress -strain- 
time  behavior  of  soil  under  various  states  of  stress  and  deformation. 

The  basic  forms  of  these  models  are  the  same  as  equation  39'.  'the  con¬ 
stitutive  coefficients,  X  ,  X  ,  p  ,  and  p^  ,  however,  are  no 
longer  constant.  One  of  the  more  general  nonlinear  versions  of  the 

Kelvin-Voigt  constitutive  equation  that  is  particularly  suited  for  soil 

9 

has  the  following  representation: 


where  oQ  ,  crB  ,  a  ,  bQ  ,  b.^  ,  a  ,  qc  ,  and  n  are  material  con¬ 

stants  that  must  be  determined  experimentally.  Equation  44  was  derived 
for  cohesionless  materials  and  includes  the  following  material 
characteristics : 

a.  Nonlinear  pressure-volume  behavior 

b.  Nonlinear  shear  stress- shear  deformation  behavior 

c.  Effect  of  pressure  on  shear  resistance 

d.  Plastic  failure 

e.  Nonlinear  viscous  effects 

A  similar  constitutive  equation  was  also  formulated  for  cohesive 
9 

materials . 
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Maxwell  model 

41.  Tne  rheological  model  of  the  Maxwell  material  consists  of  a 
linear  spring  and  a  dashpot  connected  in  series.  The  constitutive  equa¬ 
tion  of  Maxwell  material  is  given  by 


(a  J, 
'  m  1 


+  X  I. ) 
m  1' 


10 


+  23  cr.  • 
m  13 


2u  d.  . 
nm  10 


(45) 


where  a  ,  X  ,  3  ,  and  11  are  material  constants  that  must  he 

m  m  m  m 

determined  experimentally.  The  one -dimensional  representation  of  equa- 

=  ^2  = 

Accordingly,  in  view  of  equations  5>  6,  and  10,  the  one -dimensional  rep¬ 
resentation  of  equation  45  becomes 


tion  45  can  be  obtained  by  setting  d^  "  d^  =  <^3 


de 


- 

m  1 


11 


de 


m  dt 


23  a  .  + 
m  11 


2(i, 
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m  dt 


(46a) 


a  J, 

m  1 


.  d£ll 
m  dt 
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(46b) 


where,  as  before,  is  the  strain  in  the  direction  of  motion  in  uni¬ 

axial  strain  configuration  and  Og2  is  the  lateral  stress  required  to 
prevent  lateral  strain.  Equations  46a  and  46b  describe  the  motion  of  a 
Maxwell  material  in  a  uniaxial  strain  condition.  These  equations  are 
coupled  differential  equations  and  must  be  integrated  in  order  to  ob¬ 
tain  stress -strain-time  relations  for  the  material. 

Standard-linear  solid 

i(2.  The  rheological  model  of  the  standard-linear  solid  consists 
of  a  Kelvin -Voigt  element  and  a  linear  spring  connected  in  series.  The 
constitutive  equation  of  the  standard -linear  solid  may  be  expressed  as 


d.  =  (XL  +  XL  +  a  J,)  6. .  +  23  a. .  +  2p  d. .  +  2u  e. .  (47) 
13  3  1  si  s  1  xj  Ks  13  ij  ij  v  " 


where  X„  ,  X  ,  a  ,  u  ,  and  u  are  material  constants  that  must  be 
S  s  s  s  s 
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determined,  experimentally.  Equation  47  can  also  be  expressed  in  one- 
dimensional  form  similar  to  the  Kelvin-Voigt  and  Maxwell  models. 
Combined  models 

43-  Various  other  combinations  of  spring  and  dashpot  models  can 
be  constructed  by  combining,  for  example,  Kelvin-Voigt  elements  with 
Maxwell  elements  in  series.  It  should  be  pointed  out,  however,  that 
the  Kelvin-Voigt,  Maxwell,  and  standard-linear  solid  models,  or  other 
combinations  of  these  models,  are  special  forms  of  equations  9  and  38 


Hysteretic  Damping  Models 


44.  The  hysteretic  damping  models  are  generally  formulated  in 
terms  of  isotropic  incremental  elastic -plastic  constitutive  equations. 
The  nonlinear  behavior  is  incorporated  into  such  constitutive  equations 
by  making  the  elastic  moduli  functions  of  the  strain  and/or  stress  in¬ 
variants.  Hysteretic  behavior  is  taken  into  account  by  using  different 
values  of  moduli  for  loading  and  unloading.  Classical  plasticity  is 
often  incorporated  into  the  model  by  specifying  a  yield  criterion  that 
effectively  serves  to  limit  the  maximum  stress  deviators  in  the  mate¬ 
rial.  As  mentioned  previously,  a  set  of  logics  or  criteria  must  also 
be  specified  for  the  hysteretic  damping  models  to  define  various  load¬ 
ing  conditions.  Several  forms  of  hysteretic  damping  models  are  pre¬ 
sently  being  used  in  stress  analysis  studies  of  hysteretic  materials. 

As  mentioned  above,  most  of  these  models  have  evolved  along  the  lines 
of  isotropic  incremental  elastic -plastic  constitutive  equations.  A 
brief  discussion  of  the  mathematical  formulation  of  isotropic  models  is 
given  in  the  following  paragraph. 


Isotropic  incremental 
elastic -plastic  models 

45-  In  the  incremental  theory  of  plasticity,  the  incremental 

strain  tensor  de, .  is  assumed  to  be  composed  of  the  plastic  strain  incre- 

e 

ment  tensor  de. .  and  the  elastic  strain  increment  tensor  de. .  ,  i.e. 

i  n  it  ' 


de. .  =  de?.  +  de?. 
10  10  10 


(48) 
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The  basic  postulate  of  the  plasticity  theory  is  the  existence  of  a 
yield  function  or  yield  criterion  F  ,  so  that  the  material  yields  when 

F(a  )  =  0  (49) 

J.  J 

In  the  absence  of  any  further  constraint,  the  material  would  then  flow 
plastically,  undergoing  plastic  as  well  as  elastic  strains.  When  the 
stresses  are  less  than  those  required  to  satisfy  equation  49 j  the  mate¬ 
rial  will  undergo  elastic  strains  only,  i.e.,  the  material  behaves  as  a 
linear  elastic  solid.  No  stress  state  that  exceeds  the  yield  criterion 
is  allowed.  These  conditions  are  summarized  as 

d'. .  =  de?.  if  F  <  0 

ij  10 

de  .  .  =  de?.  +  de?.  if  F  =  0  ^ 

io  io  i(i 

F  >  0  hot  admissible 

46.  The  elastic  relation  of  stress  and  strain  was  stated  pre¬ 
viously  in  terms  of  the  strain  energy  functions  II  and  (equations 
16  and  20).  In  terms  of  stress  and  strain,  the  elastic  relation 
becomes 


10 


1  ±  v  v  * 

E  ij  ‘  E  °kk  ij 


from  which  one  obtains  the  elastic  strain  increment  tensor 


(51) 


1  +  v 
E 


s-  da, ,  b .  . 
E  kk  io 


(52) 


The  plastic  strain  increment  tensor  is  derived  from  the  yield  function 

10 

(equation  49)  by  utilizing  the  concept  of  the  plastic  potential,  i.e. 


10 


(53) 
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in  which  A  is  a  positive  factor  of  proportionality  that  must  be 
determined  from  the  condition  that  the  stresses  must  not  violate 
equation  49 •  Consequently,  A  is  not  unique  and  has  different  formsi 
for  various  yield  functions.  In  view  of  equations  52,  53,  and  48,  the 
total  strain  increment  tensor  takes  the  following  form 


1  +  v 
E 


E  dakk6ij 


+  A 


&F 

C )<T.  . 

U 


454) 


Equation  54,  in  conjunction  with  the  conditions  given  by  equation  5D, 
constitutes  the  incremental  stress -strain  relation  of  the  theory  of 
plasticity. 

47.  The  most  widely  used  incremental  elastic -plastic  model  an 
engineering  is  the  Prandtl-Reuss  constitutive  equation,  which  utilizes 
the  well-known  von  Mises  yield  condition  given  by 


S  S 
ran  mn 


=  k 


(55) 


where 

Smn  =  cornPonetvts  o;4  stress  deviator  tensor 
k  =  material  constant 

48.  In  view  of  equation  55?  the  yield  function  P  for  Prandtl 
Reuss  material  becomes 


F  “  I  <WU  -  *a 


(56) 


Employing  equation  56  in  equation  53,  the  plastic  strain  increment  ten¬ 
sor  becomes 


From  equation  57,  it  follows  that 


il'ii  =  ASit  =  0 


(57) 


(58) 
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indicating  that  no  plastic  volume  change  can  occur  in  the  plastic  range 
for  the  Prandtl-Reuss  material.  The  proportionality  factor  A  for  the 
Prandtl-Reuss  material  is  given  by 


where  W  is  the  rate  of  work  done  in  shear  or  the  increment  of  inter - 
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nal  energy  due  to  distortion.  Accordingly,  ir>  view  of  equations  59, 
57»  and  54,  the  strain  increment  tensor  becomes 


deij  "  1  E  V  dai$  "  E  d0kk6ij  +  A  Sij 


Equation  60  is  the  Prandtl-Reuss  constitutive  equation;  it  applies  in 
the  plastic  range  when  W  >  0  and  equation  55  holds.  If  W  <  0  ,  un¬ 
loading  is  taking  place  and  the  behavior  of  the  material  is  governed  by 
the  elastic  constitutive  equation  in  the  incremental  form,  i.e.  equa¬ 
tion  52.  If  W-'*  0  ,  the  loading  is  said  to  be  neutral,  and  equations 
52  and  60  become  identical.  The  above  results  can  be  summarized  in  the 
following  form: 


deij  =  E  V  dcij  '  E  d<7klceij 


~  S  S  <  k 
12  mn  sin 


elastic  loading 


(6la) 


_l  +  v  v  k  .  V 2  SmnSmn  ~  k 

d£i,1  =  E  d<Jij  ’  E  d0kk5i.i  *  ~2  Sij  When  $>0 


1  +  V  .  V  .  - 

de  =  — — ■  -  do.  .  -  -  da, ,  6 ,  . 
■tj  E  ij  E  kk 
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W  <  0 

when  (6lc) 

.  elastic  unloading 


From  equation  6l,  it  is  apparent  that  the  Prandtl-Reuss  material  ex¬ 
hibits  energy  dissipation  in  shear.  However ..  in  view  of  equation  58, 
the  material  behaves  elastically  during  hydrostatic  deformation. 

49 «  Other  forms  of  elastic -plastic  constitutive  equations  are 
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also  available  that  are  more  general  than  the  Prandtl-Reuss  material. 
For  example,  by  utilizing  the  Drucker -Prager  yield  condition 


in  equation  53 ?  and  utilizing  the  W  criterion,  a  more  general  consti¬ 
tutive  equation  can  be  formulated  that  can  account  for  energy  dissipa¬ 
tion  in  shear  which  is  accompanied  by  volume  expansion.  It  is  also 
possible  to  formulate  a  work -hardening  plastic  model  by  constructing  a 
yield  function  that  is  dependent  on  the  plastic  strain  tensor,  i.e. 


(63) 


The  yield  function  then  changes  or  moves  as  plastic  strain  takes  place 
in  order  to  account  for  the  work -hardening  behavior  of  the  material. 
Variable  moduli  models 

50.  The  elastic -plastic  models  fail  to  reproduce  the  stress - 
strain  behavior  of  highly  hysteretic  materials  such  as  soil  under  com¬ 
plex  and  transient  states  of  stress  and,  therefore,  are  limited  in 
their  application.  In  an  attempt  to  overcome  this  deficiency,  an  alter¬ 
nate  approach  was  developed  that  has  led  to  the  formulation  of  a  series 
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of  constitutive  models  known  as  variable  moduli  models.  Unlike  the 
elastic -plastic  models,  there  is  no  explicit  yield  function  or  flow 
rule  associated  with  the  variable  moduli  models. 

51.  The  basic  constitutive  relation  of  the  variable  moduli  model 
is  given  by 


dcr.  .  =  Kde, ,  6.  . 
10  kk  10 


2p  (  dCj } 


I  \ 

-  ~  deM  6. . 
3  kk  10/ 


where  d  =  an  incremental  change  in  stress  or  strain  ru.  1  K  =  the  bulk 
modulus  .  Equation  64  can  also  be  writcen  in  terms  of  the  deviatoric 
and  spherical  components  of  deformation,  i.e. 


1 
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deI*6i3, 


(65  a) 
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dokk  “  3Kd€kk 


(65b) 


Nonlinear  behavior  can  be  incorporated  in  the  model  by  making  K  and 
ji  functions  of  strain  and/or  stress  invariants.  The  hysteretic  be¬ 
havior  of  the  material  is  taken  inuo  account  by  using  different  values 
of  K  and  p  on  loading  and  unloading. 

52.  To  date,  in  constructing  various  variable  moduli  models,  it 
has  been  assumed  that  volumetric  strains  are  caused  only  by  changes  in 
pressure,  i.e.,  that  there  is  no  coupling  between  volumetric  strain  and 
the  components  of  the  deviatoric  stress  tensor.  Consequently,  the  bulk 
modulus  is  only  related  to  pressure  or  the  volumetric  strain  and  is 
independent  of  stress  deviations. 

53.  From  equation  64  (or  equation  65),  it  is  apparent  that  the 
constitutive  equations  of  variable  moduli  models  are  a  set  of  first- 
order  differential  equations.  The  integration  of  the  differential 
equations,  for  a  given  stress  path  and  initial  conditions,  leads  to 
stress -strain  relations.  Therefore,  a  relation  between  stress  and 
strain  is  the  outcome  of  the  theory,  but  it  depends  completely  on  the 
stress  path  used  to  reach  the  final  state  as  well  as  the  initial  state 
of  the  material. 

54.  The  most  recent  version  of  a  variable  moduli  model  was  de- 

ip 

veloped  by  Nelson.  The  load/unload  criterion  associated  with  the 
model  is  composed  of  the  criterion  for  the  deviatoric  part  of 

deformation  and  for  the  hydrostatic  part. 

55.  For  initial  loading,  the  bulk  modulus  takes  the  form 


K  -  ho  -  K0  *  %  <ekk>  +  *2  <£kk>‘ 


and  for  unloading  and  reloading 


K  =  %1/BS  "  K0u  + 


where  LD  =  loading  ,  UN  =  unloading  ,  RE  =  reloading  and  p  = 
pressure.  The  parameters  KQ  ,  ,  Kg  ,  K  ,  and  are 


31 


material  constants  that  must  be  determined  experimentally.  The 
unloading/reloading  bulk  modulus  is  used  whenever  the  pressure  is  de¬ 
creasing  or  whenever  pressure  is  increasing,  but  its  magnitude  is  less 
than  the  maximum  previous  pressure. 

56.  The  expressions  for  shear  modulus  consist  of  two  parts.  At 
low  pressure  levels,  i.e.  for  p  less  than  some  critical  pressure  pc  , 

i*ld  “  Mo  *  <Wt-  y. +  °y +  °y?2  <68a> 

MulJ  ■  Mou  +  “lu S  SmnSmn  +  °W>  +  V*  (68b) 

For  p  >  pc  ,  the  shear 

a  function  of  *  -  S  S 
M2  mn  mn 


modulus  becomes  independent  of  pressure  and  is 
only.  Accordingly,  in  view  of  equation  68 
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Pc  "  "  2c*, 
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(69) 


and  for  p  >  p  ,  the  expressions  for  shear  modulus  become 
1  c 

|  PLD  =  ^0  '  5  +  >/  2  SmnSi 
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As  indicated  by  equations  69  and  ’JO,  the  expressions  for  p 
are  continuous  at  P  “  Pc  •  The  parameters  ,  a,  , 


(70a) 

(70b) 
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UN 


"Ou  ’  '"lu  »  alu  5  and  a2u 
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ere  material  constants  that  must  be  de- 


^Ou  9  '"lu  9  "lu 
termined  experimentally.  The  loading  shear  modulus  is  used  whenever 

TQCt  >  0  during  both  initial  loading  in  shear  and  subsequent  reloading 
in  shear.  The  unloading  shear  modulus  is  used  whenever  tqc^.  <  0  . 

The  model,  therefore,  makes  no  distinction  between  initial  loading  and 
subsequent  reloading  in  shear.  This  logic,  however,  results  in  exces¬ 
sive  strains  during  cyclic  leading.  To  overcome  this  deficiency,  one 
possibility  would  be  to  use  the  unloading  shear  modulus  whenever 


t  ,  <  0  or  whenever 
oct  — 


toc^.  >  0  but  the  magnitude  of  t  ^ 


is  less 


than  its  previous  maximum  value.  This  logic  is  analogous  to  the  logic 
used  for  hulk  modulus.  Another  possibility  is  to  use  a  linear  combina¬ 
tion  of  and  for  reloading  in  shear.  The  reloading  problem, 

in  general,  is  not  well  defined,  and  additional  cyclic  data  or  theoret¬ 
ical  considerations  are  necessary  to  explain  this  phenomenon. 

57.  The  variable  moduli  models  dissipate  energy  under  both  the 
deviatoric  and  the  hydrostatic  stress  conditions  and  have  the  capabil¬ 
ity  of  matching  experimental  data  quantitatively.  There  is,  however, 

a  basic  theoretical  objection  for  their  utilization  in  the  solution  and 
formulation  of  boundary -value  problems  for  hysieretic  materials  under 
cyclic  loading  conditions.  This  objection  stems  from  the  load/unload 
criterion  for  shear  associated  with  these  models  and  the  neutral  load¬ 
ing  conditions  defined  by  roct  »  0  .  As  was  pointed  out  previously, 
whenever  tqc^  =  0  ,  the  constitutive  relations  for  initial  loading 
and  unloading  must  become  identical  in  order  to  have  unique  solutions 
for  the  boundary-value  problems  of  interest.  From  equation  68,  it  is 
obvious  that  the  expressions  for  shear  modulus  do  not  become  identical 
whenever  -  0  and  one  can  use  either  pLD  or  at  such  a 

ioad/unload  interface.  This  theoretical  objection  for  the  use  of  vari¬ 
able  moduli  models  has  led  to  the  development  of  a  series  of  elastic - 
plastic  work -hardening  constitutive  models  that  satisfy  the  continuity 

conditions  associated  with  neutral  loading  and  at  the  same  time  match 
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the  experimental  data  quantitatively. 

58.  The  constitutive  equations  of  Prandtl-Reuss  material  and  the 
variable  moduli  models  are  only  special  forms  of  the  hypoelastic  con¬ 
stitutive  equation  (equation  25). ^ 


( 


PART  V:  DISCUSSION  AND  RECOMMENDATIONS 

Discussion 

59-  The  hysteretic  character  of  the  stress-strain  behavior  of 
earth  materials  can  be  described  by  phenomenological  constitutive  <  na¬ 
tions  using  two  alternate  procedures.  In  the  first  procedure,  the  hys¬ 
teretic  behavior  of  earth  materials  is  accounted  for  by  using  different 
sets  of  constitutive  relations  for  loading  and  unloading,  in  conjunc¬ 
tion  with  a  set  of  criteria  that  define  various  loading  conditions,  i.e. 
initial  loading,  unloading,  and  reloading.  The  constitutive  equations 
in  this  case  are  time -independent.  In  the  second  procedure,  a  single 
constitutive  equation  is  specified  and  the  energy-dissipation  proper¬ 
ties  of  earth  materials  are  accounted  for  by  the  appearance  of  the  time 
derivatives  in  the  constitutive  equation,  which  in  this  case  is  time- 
de pendent. 

60.  Either  of  the  above-mentioned  procedures  can  be  used  to 
solve  soil  dynamics  boundary- value  problems  and  will  yield  stress- strain 
curves  that  exhibit  energy-absorbing  h/steretic  behavior.  The  mecha¬ 
nisms  of  energy  dissipation  in  the  two  procedures,  however,  are  differ¬ 
ent.  Additional  experimental  dynamic  data  and  theoretical  studies  in 
the  areas  of  nonlinear  viscoelasticity  and  plasticity  are  required  in 
order  to  determine  which  of  the  two  procedures  is  more  appropriate  for 
describing  the  stress -strain  behavior  of  earth  materials.  At  the  pres¬ 
ent  time,  for  reasons  of  practicability  and  depending  on  the  type  of 
problem  being  considered,  some  investigators  use  the  first  procedure 
while  others  prefer  the  second.  A  subsequent  report  (Report  2)  will 
review  the  current  laboratory  methods  of  determining  the  damping  capac¬ 
ity  of  soil;  it  will  describe  specific  techniques  used  to  evaluate 
hysteretic  effects  by  both  the  time -i idependent  and  the  time -dependent 
approaches. 


Recommendations 

6l.  Research  efforts  presently  under  way  at  the  U.  S.  Army 
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Engineer  Waterways  Experiment  Station  and  elsewhere  are  aimed  at  devel¬ 
oping  material  models  that  exhibit  energy-absorbing  hysteretic  behavior 
for  use  as  input  to  ground  motion  calculation  computer  codes.  The 
technique  used  for  formulating  the  constitutive  equations  follows  the 
procedure  for  constructing  hysteretic  damping  models.  The  constitutive 
equations  are,  therefore,  time -independent  and  do  net  manifest  viscos¬ 
ity  effects.  Very  little  effort  has  been  expended  in  recent  years  to 
formulate  mathematical  material  models  that  account  for  the  energy- 
dissipation  properties  of  the  earth  materials  through  viscosity  mecha¬ 
nisms.  These  types  of  time -dependent  material  models  are  highly  desir¬ 
able  for  the  solution  of  boundary -value  problems  involving  steady-state 
cyclic  loading  conditions.  Moreover,  they  have  the  added  advantage  of 
being  able  to  account  for  strain -rate  effects. 

62.  Thus,  it  is  recommended  that  research  efforts  to  develop 
time -dependent  constitutive  equations  for  earth  materials  for  the  for¬ 
mulation  and  solution  of  soil  dynamics  boundary- value  problems  continue. 
Such  relations  should  be  capable  of  qualitatively  and  quantitatively 
matching  the  salient  nonlinear  and  hysteretic  response  characteristics 
of  earth  materials,  not  only  as  determined  in  the  one-dimensional  con¬ 
figuration  but  also  under  a  variety  of  other  laboratory  test-boundary 
conditions.  In  the  author’s  opinion,  such  efforts  should  be  governed 
by  the  following  criteria: 

a.  The  constitutive  equations  must  take  into  account  non¬ 
linear  pres sure -volume  behavior,  nonlinear  shear  stress- 
shear  deformation  behavior,  effects  of  hydrostatic 
stress  on  shear  deformation,  shear  fracture,  spalling, 
and  plastic  failure  of  earth  materials. 

b.  The  numerical  values  of  the  various  coefficients  defin¬ 
ing  the  constitutive  models  should  have  physical  signif¬ 
icance  in  terms  of  compressibility,  shear  strength, 

etc . ,  so  that  when  extrapolating  data  for  the  different 
materials,  rational  engineering  judgments  can  be  made  as 
to  the  relative  magnitudes  of  the  constitutive  coeffi¬ 
cients  based  on  geologic  descriptions,  mechanical 
properties,  and  other  conventional  indices.  The  consti¬ 
tutive  coefficients  should  not  be  merely  a  set  of 
numbers  generated  through  a  trial -and-error,  black  box 
routine  to  fit  a  given  set  of  data. 
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63.  The  constitutive  models  must  be  used  to  obtain  analytical  or 
closed-form  solutions  to  special  problems,  e.g.  one -dimensional  wave 
propagation,  in  order  to  study  the  dissipation  characteristics  of  the 
models  as  well  as  the  theoretical  implications  of  the  new  phenomena 
they  present,  such  as  the  requirements  imposed  by  uniqueness  and  con¬ 
tinuity  considerations . 

64.  The  nonlinear  Kelvin-Voigt  model  discussed  by  Rohani  (equa¬ 
tion  44),  or  an  extension  of  this  model,  can  be  used  as  a  starting 
point  for  such  an  undertaking. 
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